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We present a detailed description of equivalent crystal theory focusing on its application 
to the study of surface structure. While the emphasis is in the structure of the algorithm 
and its computational aspects, we also present a comprehensive discussion on the calculation 
of surface energies of metallic systems with equivalent crystal theory and other approaches. 
Our results are compared to experiment and other semiempirical as well as first-principles 
calcualtions for a variety of fee and bcc metals. 
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I. INTRODUCTION 


Computational material science or the theory of materials has recently come of age. 
Calculation of properties of real materials at the atomic level such as grain boundary or dis- 
location energies or the dynamics thereof, which in the recent past have seemed intractable, 
now have some hope for realistic modelling. An even more startling assertion is that mod- 
elling of tribological phenomena is also now feasible. Tribology is a particularly difficult 
field for theoretical studies at the atomic level because of the different possible materials in 
contact under high loads, often with high degrees of disorder, but as Landman et al [1] have 
shown, a great deal of progress has been made in approaching problems of practical interest. 

More generally, there has been tremendous progress in several areas in the field of material 
science, due largely to the advent of powerful theoretical techniques and a rapid growth 
in computational capabilities. As a result, the need for a ‘theoretical laboratory’ where 
materials can be modelled, and their properties studied, at will, has become a fundamental 
part of any industrial or scientific endeavour. In order to satisfy this need, however, it 
is of essential importance that the tools needed to implement such goal have to fullfill a 
number of conditions which not necessarily complement each other. On the one hand, 
physical accuracy is of paramount importance. In this sense, first-principles approaches are 
impressively successful, in that accurate agreement with experiment was obtained in several 
jnctar w-g-g (band structure, transport and magnetic properties of solids, etc.). Properties 
of interest to the material scientist, however, are generally related to ‘defects’ in atomic 
structure which involve a partial loss of periodicity, and thus requiring calculations over a 
large number of atoms (dislocations, grain boundaries, interfaces between different materials, 
etc.). This poses an almost insurmountable obstacle, in that current technology prevents 
ab-initio calculations to have an acceptable level of computational efficiency to deal with 
problems dealing with a large number of atoms. For example, it is only recently became 
possible to predict which simple structure (fee, bcc, hep) for an elemental metal has the 
lowest energy [2]. 
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On the other hand, in order to be able to perform realistic modelling within the con- 
straints of the current available computer technology, the use of pair-potentials (two-body 
forces) offers an efficient alternative, but with a great deal of loss of physical accuracy, 
making them unreliable for predicting real material properties. 

In the last decade, due to the obvious need to fill the gap between the two alterna- 
tives mentioned above, the use of semiempirical approaches emerged as the intermediate 
choice: based on sound physical concepts, they provide generally accurate results and with 
considerable computational efficiency over first-principles approaches. In general, semiem- 
pirical approaches tackle the many-body problem by determining a functional form for the 
cohesive energy based on some physical model, including some parameters determined by 
fitting to experimental properties. Once these constants are determined, the form is used 
to calculate the energy or dynamic behavior or other properties of interest, such as defect 
energies. While sometimes oversimplifying the problem, by the global treatment that is 
made of the electronic structure of the solid, the semiempirical approaches have proven to 
be extremelly effective in accomplishing its goals. A review of recent efforts dealing with cal- 
culation of properties of materials using semiempirical techniques shows that three methods 
dominate the research activity in this area: effective medium theory [3] (EMT), embedded 
atom method [4] (EAM) and equivalent crystal theory (5](ECT). 

In this paper, we concentrate on equivalent crystal theory for the study of surface prop- 
erties. More specifically, we present a detailed application of equivalent crystal theory to the 
calculation of surface energies and surface relaxations. Some previous studies [5,6] indicate 
that ECT provides, in general, a very good description of such surface properties, based on 
comparison with available experimental data. Most impressively, ECT was recently used to 
provide a simple description of the phenomenon of avalanche in adhesion [7]: at a certain 
separation , two approaching surfaces snap together, the critical distance for such event 
depending on the thickness of the metals in contact. This example demonstrates that cal- 
culational methods such as ECT enable examination of phenomena which are difficult to 
observe experimentally. Other recent application of ECT to surface phenomena include sur- 
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face reconstruction [8], modelling of atomic force microscopy [9], friction [10] and energetics 
of steps and kinks in metallic surfaces [11]. 

We organize this paper as follows: Section 2 provides a brief summary of ECT while 
Section 3 focuses on the application of the formalism to the calculation of surface energies 
and relaxations. Section 4 reviews the results obtained for a variety of metallic surfaces. We 
devote Section 5 to the discussion of advantages, limitations and general considerations for 
future applications of ECT for the computer modelling of materials and their properties. 
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II. EQUIVALENT CRYSTAL THEORY 


Equivalent crystal theory [5] is based on an exact relationship between the total energy 
and atomic locations and applies to surfaces and defects in both simple and transition 
metals as well as in covalent solids. Lattice defects and surface energies are determined via 
perturbation theory on a fictitious, equivalent single crystal whose lattice constant is chosen 
to minimize the perturbation. Ferrante, Rose and Smith [12] found that the cohesive energy 
of metals and covalent solids as a function of varying the lattice parameter isotropically 
could be represented by a simple relationship, 

E = E c E*{r ■*) (1) 

where Ec is the cohesive energy and E*(r*) is a universal shape for the curve, which can be 
well represented by the Rydberg function 

E m (r m ) = -(1 + r*)e -r * (2) 


and 


r = 


r — rwsE 

l 


( 3 ) 


and where r is the Wigner-Seitz radius, rwsE is the corresponding equilibrium value and / 
is given by 

E c x 1/2 


1 (( <PE/dr 2 ) 


— ) 
\rwSE / 


( 4 ) 


which can also be written in terms of known experimental quantities as 

■-/ 


Ec 


( 5 ) 


12x Bt WS e ’ 

where B is the bulk modulus. The energy of the equivalent crystal, as a function of its 
lattice constant, can therefore be given by the universal binding energy relation (UBER) 
given by eq.(l). Fig. 1 displays total energy curves as a function of interatomic spacings for 
bimetallic adhesion and cohesion in metals as well as the UBER. 
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Let’s consider a certain, arbitrary defect. Let e be the total energy to form the defect or 
surface, then 

* = («) 

t 

where is the contribution from an atom i close to the defect or surface. ECT is based on 
the concept that there exists, for each atom t, a certain perfect, equivalent crystal with its 
lattice parameter fixed at a value so that the energy of atom i in the equivalent crystal is also 
d. This equivalent crystal differs from the actual ground-state crystal only in that its lattice 
constant may be different from the ground-state value. We compute £,- via perturbation 
theory, where the perturbation arises from the difference in the ion core electronic potentials 
of the actual defect solid and those of the effective bulk single crystal. 

For the sake of simplicity, the formal perturbation series is approximated by simple, 
analytic forms which contain a few parameters, which can be calculated from experimental 
results or first-principles calculations [5]. Our simplified perturbation series for £, is of the 
form 

£ , = e c \f- K(i)l + E r M'.j)] + £ *)1 + £*" \ P) 

l i i* r* J 

where 

F m [a*] = 1 - (1 + o*)e-“* (8) 

EcF m (a *) uses the UBER described above to determine the change in energy of the equiv- 
alent crystal above the ground state value. In this expression, we distinguish four different 
contributions to the energy of atom t and thus, the existence of four different equivalent 
crystals which have to be determined for each atom t. 

The first term, F* [aj(t)], contributes when average neighbor distances are altered via 
defect or surface formation. It can be thought of as representing local atom density changes. 
In most cases, this ‘volume’ term is the leading contribution to £, and in the case of isotropic 
volume deformations, it gives £,■ to the accuracy of the UBER. The value of flj(t), the 
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lattice parameter of the first equivalent crystal associated with atom i , is chosen so that 
the perturbation (the difference in potentials between the solid containing the defect and 
its bulk, ground-state equivalent crystal) vanishes. Within the framework of ECT, this 
requirement translates into the following condition from which is determined. 

NR p e~ aRl + Mi?£e- (a+ * )Rj - £ rje~ l “ +5(ri)]rj = 0 (9) 

defect 

where the sum over the defect crystal or surface is over all neighbors within second-neighbor 
(NNN) distance, r, is the actual distance between atom t and a neighbor atom j , N and 
M are the number of nearest-neighbor (NN) and next-nearest-neighbors, respectively, of the 
equivalent crystal (12 and 6 for fee, 8 and 6 for bcc) and p, a and A are parameters known 
for each atomic species. Table 1 displays the values of these parameters for the fee elements 
used in this work (see Ref. [5] for a complete list). Sfo) is a screening function and R x 
and R 2 are the NN and NNN distances in the equivalent crystal. The equivalent lattice 
parameter, a*, is thus related to the scaled quantity ^ via 

aj = — rivscj /l, (1®) 

where r WSE is the equilibrium Wigner-Seitz radius, / is a scaling length and c is the ratio 
between the equilibrium lattice constant and rwsE ■ Thus, the determination of the energy 
for an ion, in or near a defect, amounts to solving a simple transcendental equation for the 
equivalent lattice parameter. 

The higher-order terms are relevant for the case of anisotropic deformations. The linear 
independence attributed to these four terms is consistent with the limit of small pertur- 
bations which we assume for the formulation of ECT. The second term, F* [a5(*»i)]> 1S a 
two-body term which accounts for the increase in energy when NN bonds are compressed 
below their equilibrium value. This effect is also modeled with an equivalent crystal, whose 
lattice parameter is obtained by solving a perturbation equation given by 

NR\e~ aRl - Nr p 0 e-° To + A 2 r p £(#, - r 0 )c-^ R >" ro) = 0, (11) 
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where fi = 4a for the metals used in this work, and R\ is the NN distance of the equivalent 
crystal associated with the deviation of NN bond length Rj from r 0 , and r 0 is the bulk 
NN distance at whatever pressure the solid is maintained (generally, r 0 is the ground-state, 
zero-pressure value). A2 is a constant determined for each metal: for fee metals A2 is given 

by 


a e A 2 /D = l/6\/2 


( 12 ) 


and 

a e A 2 /D = l/4\/3 (13) 

for bcc metals. The quantity D is given by 

D = Ne- aro {ar 0 -p) (14) 

(see Table 1). The scaled equivalent lattice parameter is then 

a% = ~ rwss) /!• (1®) 

The third term, F m [03(1, j, k)\ accounts for the increase in energy that arises when bond 
angles deviate from their equilibrium values of the undistorted single crystal. This is a 
three-body term and the equivalent lattice parameter associated with this effect is obtained 
from the perturbation equation 

NRle~ aRl - Nrge- aro + A 3 r^e~ a{R}+Rk ~ 2ro) sm(6 jk -9) = 0 (16) 

where A3 is a constant listed in Table 1 and 6j k is the angle between the NN distances Rj 
and R k with the atom * at the center. 9 is the equilibrium angle, 70.5° for bcc and 90° for 
fee. This term contributes only when there is a bond-angle anisotropy (9 Jk ± 9). The scaled 
lattice parameter is then 

<13 = (^- — /!♦ 
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The fourth term, F* [a^ < 7 )]? describes face diagonal anisotropies (see Ref. [5] for a 
detailed description, for each lattice type, of the structural effect associated with this term). 
The perturbation equation reads 

NR\c ~ aR ' - Nr*e- Qr ° + A 4 rg l<fp = 0 (18) 

where d is the face diagonal of the undistorted cube and A 4 is a constant adjusted to 
reproduce the experimental shear elastic constants (Table 1). Finally, — t'wse'J /l- 
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III. SURFACE ENERGY CALCULATION 


A. Rigid surface 

In this section, we illustrate the formalism described in the previous section, by pro- 
viding a detailed application of ECT to the calculation of the surface energy of a metallic 
surface. Before we go through the details, we first make a few considerations of general 
interest regarding surface energy calculations and the simplifications introduced by certain 
assumptions. 

Consider a rigid surface (i.e., no interlayer relaxation): all bond lengths and angles 
retain their bulk equilibrium values, thus F*(a$) = F*(a$) = F*(al) = 0. The surface 
energy is therefore obtained by solving for the ‘volume’ term represented by F*«) only. 
We note that for this case only, a single atom need be considered per plane since each site 
is equivalent and only a few planes should be included since convergence to bulk values 
occurs rapidly. If we consider a rigid displacement of the surface layer towards the bulk, 
as is the case in most metallic surfaces, the higher-order terms become finite: some bonds 
are compressed, contributing to F*(a^), the bond angles near the surface are distorted as 
well as the difference between face diagonals in some cases, generating an increase of energy 
v j a /’•(a*) and F^aJ). For the case studied in this work, these additional contributions to 
£< are generally small, representing only 1 % to 2 % of the total energy. However, while 
these anisotropy terms are small for metals when there is no reconstruction, they play an 
important role in the energetics of these defects where the differences in energy between the 
rigid and relaxed configurations axe also small. In what follows, we will refer to this ECT 
formalism as ECT II [5]. 

As noted, the bond-angle and face-diagonal anisotropies are of very little interest for 
the calculation of metallic surface energies. For this case only, an earlier version of ECT 
[13], which we will refer to as ECT I, provides a simpler, although a bit less accurate 
framework for a surface energy calculation. The second term in eq.(7) is replaced by a 
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simple expression which allows for the direct calculation of the energy associated with bond- 
compression effects, 

N, M„ a 

t.=&EE 7=-n<cj («) 

n=l m=l ljmn 

where N, is the number of atoms in the solid, @mn — 1 if a mn — ® and ^mn — 0 otherwise, 
M n is the number of nearest-neighbors of atom n, L mn is the number of nearest-neighbors 
of atom m or n, whichever number is smaller, and a* mn is given by 


a 


Rmn/ci — rwsE 
l 


( 20 ) 


with 


1 = 


I 


Ec 

12lT BrwSE ' 


( 21 ) 


B is the bulk modulus of the crystal, Rm n is the distance between atoms m and n, ci is the 
ratio of the equilibrium nearest-neighbor distance in the crystal to rwsE , and rwsE is the 
equilibrium Wigner-Seitz radius. In ECT I [13], the third and fourth term of the energy 
expansion (eq.(6)) are ignored. In what follows, we will list results as obtained with either 
one of the two versions of ECT. Those results obtained with the full energy expansion (ECT 
II) [5], will be analyzed in terms of the value of the parameter which dictates the ‘strength’ 
of the bond compression term therefore playing an important role in the energetics of surface 
relaxation as it will be seen that this term is mainly responsible for avoiding the collapse of 
the top layers onto each other. 

We now concentrate on the calculation of rigid surface energies with ECT. Each atom 
in each plane is equivalent to all other atoms in the plane, therefore it is only necessary to 
solve the ECT equations for a single atom in each plane and the only remaining question 
is how many planes are needed for convergence. Let’s assume that the planes are labeled 
by an index j, j = 1, where j = 1 indicates the surface plane and j = b denotes a 
certain ‘bulk’ plane: a plane below which no surface effects are considered. For a single 
non-equivalent atom in plane j, we write the ECT equation 
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( 22 ) 


where N and M are the number of nearest and next-nearest neighbors, respectively, in 
the perfect crystal ( N = 12, M — 6 for fee, N = 8, M = 6 for bcc, etc.). If there is no 
change in the interatomic distance between any pair of atoms, then all nearest neighbors and 
next-nearest neighbors in the defect crystal are separated by the same distance (ri and r 2 , 
respectively), therefore, the last term in eq.(9) is greatly simplified as r, = n for the nearest 
neighbors and rj = r 2 for the next-nearest neighbors. Moreover, the screening function 
adopts the values 5(ri) = 0 and S(r 2 ) = 1/A, therefore avoiding further calculations. The 
values of Nj and Mj are easy to determine for each lattice type. Table 2 indicates some values 
for low index faces of fee and bcc metals. We stress the point that all these simplifications 
are a consequece of assuming a rigid lattice, otherwise, each interatomic distance should be 
considered separately. 

Finally, after eq.(22) is solved for all values of j, the corresponding solutions are 
used to evaluate the contribution to the surface energy a jz 

a j = E c F\a)), ( 23 ) 

where 

a » = (M j) /c~ r tvsg) (24) 

and where F* is given by eq.(8). The right hand side of eq.(22) is a measure of the defect 

(in this case, a surface). If the atom under consideration was a bulk atom in a perfect 

crystal, then Nj = N and Mj = M, which states the obvious fact that the equivalent crystal 
coincides with the real crystal. Therefore, the quantity 

Q e = Niie~ ar ' + Mr£e"<" + * )r2 , (25) 

where rj and r 2 represent the equilibrium nearest- and next-nearest-neighbor distances, 
respectively, can be considered as the reference magnitude for any defect. If we denote 
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with Q the r. h. s. of eq.(25), we then have three possibilities for any defect: a) if Q = Q e , 
as noted above, the equivalent crystal for the atom in question is the equilibrium crystal 
itself, therefore, no additional calculation is necessary and that atom does not contribute any 
additional energy to the total energy of the defect crystal. It should be noted that Q = Q e 
does not mean that there is no defect present: it essentially means that the atom in question 
finds itself in a perturbed electron density that, in average, matches the one it would find 
in a perfect crystal. Therefore, we can say that the ground state of the equivalent crystal is 
‘degenerate’, in the sense that many defects can be associated with such equivalent crystal. 
This is true, of course, for any other value of Q. b) if Q > Q e , the defect represents a 
situation where the electron density is higher than in the perfect crystal. The corresponding 
equivalent crystal must have, as a consequence, a smaller lattice parameter than the perfect 
crystal. Conversely, c) if Q < Q e then the lattice parameter of the equivalent crystal is larger 
than the equilibrium value. Examples of the last two cases are, respectively, an interstitial 
atom and a vacancy. 

As equivalent crystal theory is based on a perturbation expansion around the equilibrium 
crystal, one would expect the solutions of the ECT equation (eq.(9)) to display reasonably 
small deviations from Ri = r 0 , the equilibrium nearest neighbor separation. As a conse- 
quence, only those solutions of eq.(9) that are in the vicinity of r 0 are acceptable. Fig. 2 
shows the graphical solution of Eq.(9): the curve indicates all the possible values of the l.h.s. 
of eq. (9), while the straight line indicates the measure of the defect, Q, while the dashed 
line corresponds to the equilibrium value Q e . Two points are indicated in the horizontal 
axis: the equilibrium nearest-neighbor distance r 0 , and the acceptable solution of Eq.(??). 
Solutions to the left of the peak of the curve should therefore be ignored. 

In conclusion, the calculation of rigid surface energies is a straightforward procedure, 
which highlights the computational simplicity associated with ECT. This case is particularly 
simple because it only involves the calculation of the first term in Eq.(7) (no bonds are 
compressed and there are no bond-angle or face-diagonal anisotropies). 
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B. Relaxed surface 


Letting the planes close to the surface relax, while introducing hardly any additional 
work in the calculation of the first term in eq.(7), it ‘turns on’ the other three contributions 
to the energy. As mentioned above, the last two terms are very small and we will, for the 
moment, neglect them. For the second term - the bond-compression term - we have, as 
discussed above, two options: the ECT I formalism [13], which introduces a simple counting 
procedure where no transcendental equation has to be solved: all that is needed is how many 
bonds between the atom in question and its nearest neighbors are compressed (r < r 0 ) and 
their length. The second option, ECT II [5], while using the same ingredients than ECT 
I, introduces a transcendental equation, eq.(ll), in the same spirit as in eq.(9), where the 
anisotropy in bond length is considered as part of the perturbation. 

For small compressions, there is no quantitative difference between the two alternatives. 
For larger compressions, the ECT II approach introduces an additional, adjustable parameter 
(0) which can be fitted to some property of the crystal that would provide a good description 
of such limiting situations. The basic qualitative difference between the two approaches is the 
fact that the ECT II approach provides a continuous contribution of the bond-compression 
energy as a function of bond length, whereas the ECT I version shows a discontinuity in the 
first derivative of the energy. While this is not a major drawback in modelling the defect, it 
could be a disadvantage when ECT is used in association with statistical techniques, where 
continuity in the derivatives of the energy could be a requirement. 

Obviously, allowing the atoms in the surface region to relax, introduces the additional 
complexity of including more non-equivalent atoms to the calculation, located in deeper 
layers. As we will see in the following sections, this is unavoidable in the case of high-index 
faces, where several planes parallel to the surface plane have to be included even before 
a bulk-like environment is reached. Moreover, these cases require a careful treatment of 
relaxations perpendicular and parallel to the surface plane. 

Many-atom effects, which axe represented in ECT by the inclusion of the three-atom 
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bond-angle anisotropy and the four-atom face-diagonal distortion terms (see the third and 
fourth term in the r.h.s. of eq.(7)) axe necessary but, in the case of surface energy cal- 
culations of metals, of very little relevance. As we will see below, they introduce a small 
correction usually of the order of 1 % of the leading term in eq.(7). That is not the case 
for semiconductors, where angular anisotropies are a significant contribution to the surface 
energy. 


C. A simple example 

For the sake of clarity, in what follows we provide a simple application to an fcc(lOO) 
metallic surface, where only the surface plane is allowed to relax. All atoms in a given plane 
are identical, therefore we only need to evaluate the contribution of one single non-equivalent 
atom per plane. Moreover, only three planes have to be considered in this calculation: atoms 
in the fourth plane and below find themselves in an equilibrium, bulk-like environment. 

Fig. 3 shows a side view of such surface: the dots indicate atoms on the plane of the 
paper and the crosses indicate atoms one plane below or above. The top plane (j = 1) is 
slightly contracted, by a distance x toward the first plane below the surface ( j = 2), but both 
the j = 2 and j = 3 planes, as well as any other plane below, are at equilibrium positions. 
The rigid interplanar spacing is d. — o e /2 (a e is the equilibrium lattice parameter), therefore, 
the distance between j - 1 and j = 2 is a e /2 - x. Following the labeling of Fig. 3 we write 
the ECT equations for the atoms A\, Bi and C u in terms of the distances r mn between 
atoms m and n: 

12 R\t- aR ' + 6i^e- (a+ * )fij = 

4r5 +4r$ lSl e-"- , '>> + 

0 = 1 ) (26) 

12 R\t~ aRl + 6B%e~ {a+ * )R ' = 

+ ^ BlAl e~ arB ^ + 5r| lB3 e-<“ + *>^ (j = 2) (27) 
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\2R\e~ aR ' + 6i^c- (a+ i )ila = 


12r p Ci c 2 e- aTc ^ + 5r^ lC3 e-(° + *> rc ^3 + r* ClAl e (a+i)rc >^ (j = 3) 


where 



rA,A 3 = r BiBi = rciC 3 =$0-a e 


and 


r Al c, = rc,^! =s x = a e -x. 

Eqs.(26)-(28) are then solved, for each value of x, for the lattice parameter 
corresponding equivalent crystal: 

a ( i\x) = 4 J) = V2 R? j = 1, 2, 3. 

Finally, the ‘volume 1 contribution to the surface energy is 

„, = £>«=. Bct,F-(c- J (x)) 
i= i i=i 

where 

- 'X 1 r i 

a U = j(— “ rwssh 

and F * is given by eq.(8). 

We now compute the bond length anisotropy contribution to the surface energy, 
the ECT I framework, the contribution is given by eq.(19): 

CT ' = E c (\r(a\ llh ) + (<■**)) 


(28) 

(29) 

(30) 

(31) 

(32) 
of the 

(33) 

(34) 

(35) 
Within 

(36) 
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where 


( 37 ) 


. . r x /ci - rwsE 

a AjBi ~ a BiAi ~ i 

as long as x > 0. If x < 0 then <rf CTI = 0. 

In order to compute the bond-compression contribution within the ECT II framework, 
we need to solve the following transcendental equation in order to find the equivalent crystals 
(for atoms Ai and Bi) associated with this defect: 

12 R\e~ aR ' - 12 rge _0,r ° + 4 A 2 Ti(r x - r 0 )e _/5(rx_,o) = 0 ( atom Ai) (38) 

The equation for atom Bu for this particular example, is identical to eq.(38) (reflecting the 
identity a BlM = a\ xBx for the ECT I case). Eq.(38) is then solved with respect to R\ and 
the energy contribution is then 

<rf CT " = E c F'(a' 2 (z)) (39) 

where a$ is given by eq.(15). 

The third term in eq.(7) deals with bond-angle anisotropies. Contributions to <r 3 come 
from atoms for which the angle between its nearest neighbors departs from its equilibrium 
value 0 O (0o = 90° for fee). Following the convention described in ref. [5] that says that if an 
atom is missing one or more nearest-neighbors, then k)) — 0 for that atom, then 

the atoms in the top layer do not contribute to cr 3 . The only contribution then arises from 
the atoms in the first layer below the surface (j = 2), for which we solve the transcendental 
equation (see eq.(16)) 

12 - 12rge- aro + A 3 rge-° (T *- ro) sin 1 0 - 0 O | = 0 (40) 


where 


*“*“ { 2r.ro "T 

Eq.(40) is solved with respect to R t and the bond-angle contribution is then 


( 41 ) 
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of CT " = iE c F-(al(x)), 


(42) 


where a| is given by eq.(17). 

For the particular case of the fee (100) face, the fourth term vanishes, as the diagonals 
associated with each atom in the layer j = 2 do not change relative to each other, although 
the length of each diagonal is different from the equilibrium value. 

Summarizing, the surface energy of an fcc(100) face, where only the top layer is allowed 

to relax is given by 

c- ECT ' = ME n-wW) + («) 

i= i 

where <jJj(cc) and by eqs.(35) and (37), respectively. Within the EOT II 

formalism, a is given by 

0 ectii = Ec{ j- F - {a - j(x) ) + F*(a;{x)) + 4F*(«;(x))} (44) 

i=i 

where ^d 03(2) 3X6 S^ ven by eqs.(15) and (17), respectively. 
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IV. RESULTS 


ECT, in any of its two versions, has been applied with a great deal of success to sev- 
eral problems involving surface phenomena [5-10]. In the case of surface energies, ECT is 
particularly successful as it gives a very good approximation to the currently accepted exper- 
imental and first-principles values, therefore making it a good candidate for predicting such 
quantity in other cases. We have carried out an extensive calculation of surface energies of 
fully relaxed (but unreconstructed) surfaces of several fee and bcc metals, and for a variety 
of choices for the Miller indices, ranging from (111) to (771) [6]. Although some of these 
surfaces are known to reconstruct, we did not look for this effect in our calculations with 
the sole purpose of generating a large data set of surface energy values from which patterns 
can be extracted, by studying the dependence of surface energy on several variables (Miller 
indeces, cohesive energy, roughness [14], borocity [15], lattice structure, etc.). 

In the previous section, we repeatedly mentioned the fact that the many-body 
anisotropies (bond-angle, face-diagonal) that appear in the ECT II formalism are not of 
great relevance for the calculation of surface energies. This point is proven by the results 
displayed in Table 3 where we show some results for low- and high-index faces, both rigid 
and fully relaxed, for the surface energy and the different contributions arising from the 
ECT II expansion (eq.(7)). Table 4 proves another point made in the previous section: in 
general, both ECT I and ECT II give results of comparable quality in the case of surface 

energies. 

Tables 5 and 6 display an extensive set of results for fee and bcc metals, where we 
compare with several experimental estimates as well as the surface energies obtained by 
other approaches. First, we examine the quality of our results by comparison with first- 
principles calculations. First-principles values do not suffer from the ambiguities in surface 
geometry and cleanliness as the difficult to obtain experimental values. The surface energy 
for polycrystalline solids is usually assumed to be close to that of the most dense (less rough) 
surface, therefore we discuss results for the densest packed fee (111) planes. We extract our 
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values from Tables 5 and 6 , for direct comparison with other methods for which we could 
find surface energies , where the first number will be the ECT II value, (), the second EAM 
[4,19], 0> the third Finnis-Sinclair [16], [Q], and the last first-principles [17], {} (all values 
axe in ergs/cm 2 ): Cu (1758), [1170,1184], {2100}; Ag ( 1221 ), [620,655], [[620]], {1221}; Pd 
(1696), [1220,1050], {1640}. We can see that in all cases ECT II agrees extremelly well with 
first-principles results, whereas Finnis-Sinclair and EAM results can be as much as 50 % 
low. 

The next comparison made is with the planar anisotropy of the surface energy as ob- 
tained by different approaches. There axe arguments [17,18] for expecting the variation in 
the surface energy from plane to plane to be smaller than expected from variations in surface 
roughness. For example, as mentioned earlier, electron smoothing in order to reduce the ki- 
netic energy could reduce the variation in surface energy expected from roughness variations. 
For this discussion, we reproduce results presented by Methfessel et al [17] (table 7) with 
additional analysis, updated values for ECT from the present paper, and more recent EAM 
results [19]. First we note that when we examine the anisotropy ratios 7110-111 = 0110/0111 
and 7 , 00-111 = 0100 / 0111 , where <r hkl is the surface energy of the face (hkl), as well as 
7100-110 = 0100 / 0110 , Methfessel et al’s assertion that the other semi-empirical methods, 
EAM and Finnis-Sinclair, show a smaller surface anisotropy than ECT is not substantiated. 
The surface anisotropies predicted by the three methods are actually quite close and all 
three predict a small variation between the values of 7100 - 110 * 

A simple modification of ECT II can bring the surface anisotropy much closer to the 
Methfessel et al values [17]. As we stated earlier, ECT introduces a screening function for 
next and higher order neighbors [5] based on the bulk screening length. The functional 
form for this screening is assumed to be an exponential which although reasonable, is ad 
hoc. Assuming no screening (Table 8 ) for both Pd and Ag, the surface anisotropy variation 
can be brought much more in line with the values reported in ref. [17], without affecting 
the outstanding agreement for the surface energies for the (111) surfaces. We note that 
such an assignment of screening length is directly comparable to that used in EAM. We 
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further note that there is not at present a consistent comparison between first-principles 
and semiempirical methods. Limiting our example to ECT, the primary input parameters 
for the method are the experimental values for the cohesive energy, the zero-pressure bulk 
modulus, the equilibrium lattice parameter and the vacancy formation energy. In order 
to make a consistent comparison it would be necessary to use input values from the first- 
principles calculations for determination of the fitting parameters in each method. 

Next we address the bcc structure metals. Our values for the ratios 7100 -no = o'ioo/^iio 
and 7111-110 = oiii/ono for W and Fe are, respectively, 7100-110 = 1-78, 1.97 and 7m_uo — 
1.32, 1.40. The values of Methfessel et al [17] and Finnis-Sinclair [16] for 7100-110 are respec- 
tively, 1.21 and 1.17 for Nb and 1.13 and 1.15 for Mo. Here, the assertion in ref. [17] is more 
accurate, however we note that for the only first-principles comparisons that we have found 
for a bcc system, that is, W, Fe (100), the ECT value of 5808 ergs/cm 2 as compared to 
5100 ergs/cm 2 for W (100) [20] (an overestimate of 14 %) and 3353 ergs/cm 2 as compared 
to 3100 ergs/cm 2 for Fe (100) [21] (an overestimate of 8 %) are in quite good agreement. 
Although the surface anisotropy seems high in this case, there is as yet no complete test of 
the bcc systems which warrant a definite conclusion. As mentioned above, we have shown 
that a simple modification of ECT can be made which corrects for the seeming discrepancy 
in surface anisotropy while maintaining the simplicity of the method and the highly accurate 
predictions for the surface energies of the densest packed planes. 

Simple bond-cutting arguments have been used to explain the face dependence of the 
surface energy, which could be thought of as the energy cost for cutting the bonds to one 
nearest-neighbor times the number of removed NN atoms in the process of surface formation. 
ECT provides a simple qualitative and quantitative description of this process: if we consider 
a rigid surface (so that second-order effects due to anisotropies are not included), one would 
reasonably expect the atoms in the top layers to be mostly affected by the creation of the 
surface. For example, for fee (100) planes, the contribution of an atom in the surface plane 
accounts for most of the surface energy, while atoms in the second layer contribute less than 
1 %. The ECT equation for determining the lattice parameter of the equivalent crystal 


21 



associated with an atom on the surface is (see eq.(22)) 

NI%e - aR 1 + = Nrfe-”' + M 1 ^e" ( “ + ^ ra , (45) 

where N\ and M\ represent the number of NN and NNN of the atom on the surface (thus, 
N-Ni and M- M x are the number of NN and NNN bonds broken in the process of creating 
the surface). Ri and # 2 axe the NN and NNN distances in the equivalent crystal of lattice 
parameter cleq- This is a transcendental equation for oeq and although it cannot be solved 
analytically, it can be shown that a simple linear approximation gives a result for a E Q which 
is only 10 % off the exact value. This approximation entails a linear dependence on the 
number of broken bonds, N - N x and M - M x . The dependence of the surface energy on 
these quantities is given by the Rydberg function, thus establishing a relationship between 
surface energy, cohesive energy and the number of NN and NNN broken bonds, 

*„rf = E c {l-(l+ «>-•}. ( 46 ) 

where a* is the scaled lattice paramater of the equivalent crystal: 

«• = q ( a *e~Jd. (47) 

The previous analysis provides a simple explanation for the differences in surface energy 
between different planes. Nevertheless, the surface energies obtained with ECT compare very 
well with experiment in all cases, as shown in Tables 5 and 6 where the ECT predictions, 
as well as other theoretical and experimental figures are displayed for fee (Table 5) and bcc 
(Table 6) metals. 

We once again note that the surfaces energies displayed are for relaxed but unrecon- 
structed surfaces. This is obviously unrealistic for high-index faces, or even cases like Au 
(100), but it provides a rough idea of the dependence of the surface energy on the roughness 
of the surface. A comparison of the results in Table 5 (obtained with ECT II allowing only 
for planar relaxations in a direction perpendicular to the surface plane) and the surface 
energies quoted in Table 4 (which also include atomic relaxations parallel to the surface 
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plane) shows that adding this degree of freedom in the calculation, while introducing a sig- 
nificant change in terms of the number of variables with respect to which the energy has to 
be minized, introduces little or no change in the values of the surface energy. 
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V. CONCLUSIONS 


Equivalent crystal theory, like other semiempirical techniques, is supposed to narrow the 
gap between formal first-principles calculations and empirical approaches. To achieve this 
goal, ECT provides a simple operational procedure with minimal computational complexity. 
At the same time, it has a strong foundation on perturbation theory. In the case of surface 
phenomena, ECT is an excellent tool for describing the different aspects of surface structure 
and, as we intended to show in this work, allows for simple visualization of the relevant 
features. 

It is an encouraging fact that, in spite of the reasonable success obtained in the ap- 
plications studied since the inception of this method, it is amenable to extensions in a 
straightforward fashion: a quick analysis of the perturbation series expansion provides clear 
indication of how the method can be perfected to deal with more complex situations. This 
flexibility is a necessary condition for any method that attempts to describe, in one single 
scheme, all possible systems. 


24 


ACKNOWLEDGMENTS 


Fruitful discussions with Dr. N. Bozzolo are gratefully acknowledged. This work was 
partially supported by the Engineering Directorate, NASA Lewis Research Center. 


25 



REFERENCES 


[1] U. Landman, W. D. Luedtke and E. M. Ringer in Fundamentals of Friction, Eds. I. 
L. Singer and H. M. Pollock, NATO ASI Series, Kluwer Academic Publishers, The 
Netherlands, 1992. 

[2] H. L. Skriver, Phys. Rev. Lett. 49, 1968 (1982); M. T. Yin and M. L. Cohen, Phys. 
Rev. B 26, (1982). 

[3] J. K. Norskov, Phys. Rev. B 20, 446 (1979). 

[4] S. M. Foiles, M. I. Baskes and M. S. Daw, Phys. Rev. B 33, 7893 (1986). 

[5] J. R. Smith, T. Perry, A. Banerjea, J. Ferrante and G. Bozzolo, Phys. Rev. B 44, 6444 
(1991). 

[6] A. M. Rodriguez, G. Bozzolo and J. Ferrante, Surf. Sci. 289, 100 (1993) and refrences 
therein. 

[7] J. R. Smith, G. Bozzolo, A. Banerjea and J. Ferrante, Phys. Rev. Lett. 63, 1269 (1989). 

[8] B. Good and A. Banerjea, Mat. Res. Soc. Symp. Proc. 278, 211 (1992). 

[9] A. Banerjea, J. Ferrante and J. R. Smith, J. Phys.: Condens. Matter 2, 8841 (1990). 

[10] G. Bozzolo, J. Ferrante and J. R. Smith, Scripta Met. Mater. 25, 1927 (1991). 

[11] S. Khare and T. Einstein (private communication) 

[12] J. H. Rose, J. R. Smith and J. Ferrante, Phys. Rev. B 28, 1835 (1983); J. H. Rose, J. 
R. Smith, F. Guinea and J. Ferrante, Phys. Rev. B 29, 2963 (1984). 

[13] J. R. Smith and A. Banerjea, Phys. Rev. Lett, bf 59, 2451 (1987); ibid., Phys. Rev. B 
37, 10411 (1988). 

[14] J. Sokolov, F. Jona and P. M. Marcus, Sol. State Comm. 49, 307 (1984). 

[15] G. Bozzolo, A. M. Rodriguez and J. Ferrante, Surf. Sci. (in press). 


26 



[16] M. W. Finnis and J. E. Sinclair, Philos. Mag. A 50, 45 (1984). 

[17] M. Methfessel, D. Hennig and M. Scheffler, Phys. Rev. B 46, 4816 (1992). 

[18] J. P. Perdew, Y. Wang and E. Engel, Phys. Rev. Lett. 66, 508 (1991). 

[19] T. Ning, Q. Yu and Y. Ye, Surf. Sci. 206, L857 (1988); T. Ning and Y. Ye (private 
communication). 

[20] C. L. Fu, S. Ohnishi, H. F. Jansen and A. J. Freeman, Phys. Rev. B 31, 1168 (1985). 

[21] J. G. Gay, J. R. Smith, R. Richter, F. J. Arlinghaus and R. H. Wagoner, J. Vac. Sci. 
Technol. A 2, 931 (1984); J. R. Smith, J. Ferrante, P. Vinet, J. G. Gay, R. Richter 
and J. H. Rose, in : Chemistry and Physics of Fracture, Eds. R. H. Jones and R. M. 
Latanision (Nijhoff, Hingham, MA, 1987). 

[22] K. -M. Ho and K. P. Bohnen, Phys. Rev. B 32, 3446 (1985). 

[23] J. A. Appelbaum and D. R. Hamann, Sol. State Comm. 27, 881 (1978). 

[24] W. R. Tyson, J. Appl. Phys. 47, 459 (1976), and references therein. 

[25] W. R. Tyson and W. A. Miller, Surf. Sci. 62, 267 (1977). 

[26] D. Wolf, Surf. Sci. 226, 389 (1990). 

[27] S. P. Chen and A. F. Voter, Surf. Sci. 244, L107 (1991). 

[28] M. Weinert, R. E. Watson, J. W. Davenport and G. W. Fernando, Phys. Rev. B 39, 
12585 (1989). 

[29] G. J. Ackland and M. W. Finnis, Philos. Mag. A54, 301 (1986); G. Ackland, G. Tichy, 
V. Vitek and M. W. Finnis, ibid. A56, 735 (1987). 

[30] S. P. Chen, Surf. Sci. 274, L619 (1992). 

[31] L. E. Murr, in Interfacial Phenomena in Metals and Alloys, Addison- Wesley, Reading, 


27 



MA, 1975. 


[32] W. R. Tyson, R. A. Ayres and D. F. Stein, Acta Met. 21, 621 (1973). 

[33] J. J. Gilman, in Fracture, edited by B. L. Averbach, D. K. Felbeck, G. T. Hahn and D. 
A. Thomas, John Wiley, NY, 1959. 

[34] A. G. Gvozdev and L. I. Gvozdeva, Fiz. Met. Metalloved 31, 640 (1971). 

[35] M. H. Richman, Trans. ASM, 60, 719 (1967). 

[36] E. N. Hodkin, M. G. Nicholas and D. M. Poole, J. Less-Common Met. 20, 93 (1970). 

[37] G. Treglia, M. C. Desjonquere and D. Spanjaard, J. Phys. C: Condens. Matter 16, 2407 
(1983). 

[38] C. T. Chan and S. G. Louie, Phys. Rev. B 33, 2861 (1986). 


28 



Element 

P 

1 

a 

A 

W*M[D 

W^aJd 

10 ~*D 

AE 


A1 


0.336 

2.105 

0.944 

7.822 

2.104 

591.4 

3.34 

4.05 

Cu 


0.272 

2.935 

0.765 

5.784 

2.530 

99.74 

3.50 

3.615 

Ag 


0.269 

3.337 

0.756 

5.390 

2.285 

12.90 

2.96 

4.086 

Au 

10 

0.236 

4.339 

0.663 

4.047 

1.673 

1.127 

3.78 

4.078 

Ni 

D 

0.270 

3.015 

0.759 

7.382 

2.793 

100.1 

4.435 

3.524 

Pd 

H 

0.237 

3.612 

0.666 

5.242 

2.012 

11.25 

3.94 

3.89 

Pt 

m 

0.237 

4.535 

0.666 

5.789 

1.727 

1.071 

5.85 

3.92 

Fe 


0.277 

3.124 

0.770 

9.183 

1.887 

60.62 

4.29 

2.86 

W 

10 

0.274 

4.232 

0.770 

12.03 

1.497 

1.179 

8.66 

3.16 


Table 1 : Computed constants and experimental input for ECT II. The constant p is 2 n — 2 , 
where n is the atomic principal quantum number, l (in A) is a screening parameter (see 
text). The constants A 3 and A 4 are dimensionless. E c (in eV) is the cohesive energy and 
a e (in A) is the equilibrium lattice constant. 


j 

(hkO 

fee 

bcc 

Ni 

Mi 

Ni 

Mi 


(100) 

8 

5 

4 

5 

1 

(110) 

7 

4 

6 

4 


(111) 

9 

3 

4 

3 


(100) 

12 

5 

8 

5 

2 

(110) 

11 

4 

8 

6 


(111) 

12 

6 

7 

3 


(100) 

12 

6 

8 

6 

3 

(110) 

12 

6 

8 

6 


(111) 

12 

6 

7 

6 


Table 2: N um ber of nearest- neighbors (Nj) and next-nearest-neighbors (Mj) of an atom in 
plane j (j = 1 is the surface plane) for (hkl) fee and bcc surfaces. 
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hkl 

Rigid 

Relaxed 

<*\ 

a 2 

<r 3 

*4 ! 


111 

1813.90 

1757.63 

1735.83 

20.48 

0.00 

1.32 


100 

2362.26 

2301.70 

2276.59 

22.61 

2.51 

0.00 


110 

2443.38 

2365.87 

2335.63 

28.32 

0.00 

1.92 

Cu 

311 

2432.55 

2363.21 

2335.13 

25.72 

1.06 

1.30 


331 

2306.33 

2237.46 

2210.38 

24.88 

0.04 

2.16 


210 

2617.94 

2543.69 

2514.30 

27.31 

0.32 

1.75 


211 

2288.15 

2227.98 

2203.65 

22.03 

0.57 

1.72 


310 

2598.22 

2528.68 

2500.59 

26.06 

0.90 

1.13 



1752.92 

1695.55 

1673.79 

20.98 

0.00 




2303.30 

2240.53 

2215.23 

23.26 

2.04 

1 



2382.47 

2303.44 

2273.40 

28.87 

0.00 

1.17 

Pd 


2368.39 

2297.26 

2269.23 

26.29 

0.93 




2242.51 

2172.17 

2145.44 

25.34 

0.03 

1.35 



2553.10 

2477.53 

2448.79 

27.41 

0.26 

1.07 


211 

2223.16 

2157.22 

2131.50 

24.02 

0.58 

1.12 


310 

2533.69 

2462.54 

2434.30 

26.79 

0.70 

0.75 


110 

1778.77 

1703.92 

1671.38 

28.43 

2.88 

1.21 


100 

3447.37 

3353.12 

3307.87 

39.89 

5.36 

0.00 


211 

2228.52 

2131.37 

2091.82 

36.28 

0.32 

2.95 

Fe 

310 

2978.04 

2866.88 

2798.64 

57.92 

8.10 

2.22 


111 

2477.70 

2390.38 

2355.99 

31.82 

0.01 

2.56 


210 

2668.29 

2583.46 

2540.82 

36.62 

4.31 

1.71 


311 

2685.95 

2578.01 

2529.75 

43.93 

2.16 

2.17 


331 

2139.92 

2050.28 

2005.81 

37.39 

3.72 

3.35 


Table 3 : Rigid and relaxed surface energies (in erg/ cm 2 ) for Cu, Pd, and Fe surfaces. Oj 
(j = 1,...,4) denotes the contribution to the surface energy from the j term in the ECT 
expansion. 


Surface 

Rigid 

ECT I 

ECT II 

fjAl 

& 0) 

1493 

1424 

1404 

/t N* 

**(210) 

3407 

3330 

3306 

Jen 

a (210) 

2618 

2559 

2540 


Table 4: Surface energies (in erg/cm 2 ) of the rigid and relaxed (210) surfaces of Al, Ni and 
Cu computed with ECT I and II. 
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Table 5 : Experimental (Exp.) and theoretical surface energies for fee metals. 
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Table6 : Experimental (Exp.) and theoretical surface energies for bcc metals. 
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■ 
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0110 
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7100-111 

1.13 



1.12 

1.19 

1.32 

Pd 

7110-111 

1.20 



1.22 

1.30 

1.36 


7ioo-iio 

0.94 

0.92 


0.92 

0.91 

0.97 


0iii 


wmm 

mm 

0.62 

0.65 

1.22 

Ag 

0100 



BBS 

0.71 

0.79 

1.58 


0110 

1.26 

Bfl 

Bin 

0.77 

0.89 

1.62 


7100 -in 

1.00 

■ 


1.14 

1.21 

1.29 

Ag 

7110 -in 

1.04 



1.24 

1.37 

1.33 


7ioo-iio 

0.96 

in 

0.94 

0.92 

0.89 

0.97 


Table 7: Comparison of ECT surface energies (in erg/cm 2 ) for Pd and Ag with results of 
full-potential LAPW calculations (ref. 17), Finnis-Sinclair (ref. 29) and EAM (refs. 4, 19) 
semiempirical results, aw denotes the surface energy of the (hkl) surface and 7 w-hkl' 
denotes the ratio between aw and a hkl' ■ 




mm 

m&m 

Ref. 17 


(100) 

1613 

1220 

■SJTS 

Ag 

(110) 

■ 

1440 
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(111) 
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(100) 

2300 

MUTIJII 

I860 

Pd 

(110) 

2380 

1 

1970 


(111) 

1750 

IkSM 

1640 


Table 8: Comparison of the surface energies (in erg/cm 2 ) of the low-index faces of Ag and 
Pd obtained with ECT and first-principles calculations (ref. 17). 
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Scaled lattice constant, a* 


Fig. 1: (a) Scaled binding energy per atom of a crystal as a function of interatomic 
separation for representative solids, (b) Scaled binding energy versus scaled separation for 


representative cases of cohesion, bimetallic adhesion, chemisorption and diatomic molecule. 
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Fig. 2: Graphical solution of eq. (9) (see text). 
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Fig. 3: Side view of a fee (100) surface. Circles denote atoms on the plane of the page, 
crosses indicate atoms one plane below or above. The index j labels the planes starting with 
j = 1 for the surface plane. 
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